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ABSTRACT 


Calculations for the compressible, equilibrium 
dissociating stagnation-point boundary layer with suc- 
tion and injection have been performed using a digital 
computer. Results are shown which reproduce both 
the perfect gas, linear viscosity -temperature re- 
lation calculations of Cohen and Reshotko, and the 
dissociating-air results of Fay and Riddell, Calcula- 
tions for the dissociated-air boundary layer with 
air injection indicate significantly higher heating 
rates than those obtained for the perfect gas, linear 
viscosity-temperature relation model. A procedure 
is indicated which allows application of the results 
for air injection to engineering calculations for 
cases of foreign-gas injection. 
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THE STAGNATION-POINT BOUNDARY LAYER WITH SUCTION 


AND INJECTION IN EQUILIBRIUM DISSOCIATING AIR 

By Kenneth C. Weston 
Manned Spacecraft Center 

SUMMARY 


Calculations for the compressible, equilibrium dissociating stagnation-point 
boundary layer with suction and injection have been performed using a digital com- 
puter. Results are shown which reproduce both the perfect gas, linear viscosity- 
temperature relation calculations of Cohen and Reshotko, and the dissociating air 
results of Fay and Riddell. Calculations for the dissociated-air boundary layer with 
air injection indicate significantly higher heating rates than those obtained for the per- 
fect gas, linear viscosity-temperature relation model. A procedure is indicated which 
allows application of the results for air injection to engineering calculations for cases 
of foreign-gas injection. 


INTRODUCTION 


Transpiration cooling has been long considered as a means for cooling the sur- 
faces of high-speed flight vehicles. Practical engineering considerations, however, 
have tended to limit its application. The emergence of the ablation heat shield for the 
protection of nose cones and manned spacecraft during reentry has provided a passive 
means of transpiring gases into the surrounding flow to alleviate aerodynamic heating. 
Thus, the calculation of the effects of gas injection into the boundary layer is now a 
problem of practical interest to the spacecraft designer. 

At orbital and superorbital flight speeds, the kinetic energy of reentry vehicles 
is sufficient to produce marked differences between the surrounding air- species com- 
position and the free- stream composition. The theory of stagnation-point heat transfer 
was extended in reference 1 to account for such composition differences resulting from 
dissociation of air molecules. It was shown there that significant quantities of energy 
may be transported through the boundary layer by the diffusion of atoms which release 
their dissociation energy upon recombination. Clearly, such a mechanism should be 
considered when sizable concentrations of atoms are present in the boundary layer. 
Concentrations of atoms and molecules on a stagnation streamline are readily esti- 
mated by solution of the normal- shock equations together with equilibrium air proper- 
ties (refs. 2 and 3). At still higher flight speeds, appreciable ionization occurs in the 
stagnation region. Such effects are not considered here. 



Calculations of the effects of mass injection into an undissociated-air stagnation- 
point boundary layer have been made by Reshotko and Cohen in reference 4. They 
found, as have other investigators, that large decreases in heat transferred to a sur- 
face may be obtained by mass transport through the surface. The present work con- 
siders the effects of suction and injection on the dissociating-equilibrium-air boundary 
layer. In essence, the equations of Fay and Riddell (ref. 1) are solved subject to the 
boundary conditions relevant to mass transfer at the wall. An engineering method of 
adapting the present results for air to the prediction of stagnation-point heat transfer 
with foreign-gas injection is also given. 

This document replaces the original NASA TN D-3889 issued in March 1967. As 
a result of an erroneous input to the computer program, the solutions, discussed under 
the heading " Equilibrium Real- Air Solutions” and tabulated in table I of the former re- 
port, were in error. The present document provides the corrected solutions as well 
as some additional results. The computer program as listed in the appendix, as well 
as the comparisons with the work of Fay and Riddell and of Cohen and Reshotko, re- 
main unchanged. 



SYMBOLS 

(pv) w 

mass-transfer coefficient, — ^ — 

'W'st 

species mass fraction 
specific heat at constant pressure 

mean specific heat, SC.c. . 

i i P> i 

ratio of specific heat of injected specie to that of air 

diffusion coefficient 

function defined in equation (12) 

function defined in equation (7) 

dimensionless enthalpy defined in equation (8) 

enthalpy 

dissociation energy 
thermal conductivity 
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pc p D 

Lewis number, — g— 

density-viscosity ratio defined in equation (13) 


Stanton number in the absence of wall mass transfer, 


Nusselt number, 


k ( h S- h w) 


static pressure 


heat -transfer rate 


heat-transfer rate at zero-mass transfer 


Reynolds number, 


cylindrical coordinate of surface normal to body axis 
dimensionless mass fraction defined in equation (9) 
temperature 

velocity parallel to body surface 
free-stream velocity 

velocity normal to body surface 
coordinate measured along body surface 
coordinate normal to body surface 


a u 2 > a 3 coefficients defined in equation (18) 


r l’ y 2 


coefficients defined in equation (19) 


similarity variable defined by equation (6) 
viscosity 

similarity variable defined in equation (5) 
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p 


density 


n 

cr Prandtl number, -^ E - 

co^ rate of generation of specie i 


Subscripts: 

e 

i 

P 

S 

w 

X 

y 

v 

$ 


denoting boundary-layer edge conditions 
pertaining to the ith specie 
denotes constant pressure 
denotes stagnation point in external flow 
denotes wall conditions 

partial derivative with respect to x coordinate 
partial derivative with respect to y coordinate 
partial derivative with respect to 77 coordinate 
partial derivative with respect to £ coordinate 


ANALYSIS 


Boundary- Layer Equations 

The present analysis employs the axisymmetric laminar compressible boundary- 
layer equations in the form derived in reference 1 for a mixture of thermally perfect 
gases. The boundary-layer equations for this case are 


(pru)„ + (prv) = 0 
A y 


a) 


puC i, x + P vC i,y 


■ CV^ut), 


+ or. 
1 


i = 1, 2, 


n 


( 2 ) 


puu x + pvu y = -P x + ( UU y) y 


( 3 ) 


4 



( 4 ) 



Subscripts x and y here indicate partial differentiation with respect to x and y. 
Terms given in reference 1 for thermal diffusion in the boundary -layer equations and 
subsequently neglected have been omitted from equations (2) and (4) since effects of 
thermal diffusion are negligible for stagnation temperatures less than 10 000° K for 
an equilibrium boundary layer. 

In the case of equilibrium at a given flight condition, thermodynamic properties 
and species-mass-fraction distributions through the boundary layer may be considered 
to be functions of enthalpy only for a specified boundary-layer-edge pressure. The 
species-conservation equation (2) is then not required and it is necessary only to work 
with equations (1), (3), and (4). 

These equations are transformed using the (£, ??) similarity variables 


{ = i We'** 


and 


n 



( 6 ) 


In addition a stream function f , dimensionless enthalpy g, and mass fraction are 
defined where 


u _ 9f 
u e = 3 V 


«n 


g = 


*4 

h s 


( 8 ) 


5 



and 


S. 


1 



(9) 


The transformed boundary- layer equations for the stagnation point become 



( 10 ) 


|-(l+d)g^ + fg 77 = 0 


(ID 


where 


d = 




(L- 1) 



P 


( 12 ) 


and 


i = 


PI-+ 

P M- 

ww 


(13) 


Equations (10) and (11) may be considered as simultaneous nonlinear ordinary differ- 
ential equations for f and g in terms of a single independent variable rj if the func- 
Po 

tions i, a, — , and d are constants or functions of 7j only. In the present work, 

P 

the Prandtl number a is taken to be a constant and the Lewis number L is taken to 

P s 

be unity. The functions £ and — are correlated as functions of g (rj). These corre- 
lations are described in a later section. 
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Boundary Conditions 


The boundary conditions relevant to the solution of equations (10) and (11) with 
suction and injection are 


f (o) = f w y°°) = i g(°) = g w 

f (0) = 0 g(°°) = 1 

The rate of mass transfer at the wall is related to f . The desired relation may 

w 

be obtained by considering the transformed continuity equation given in reference 1 as 


pv= 




(14) 


VF f. 

Noting that — p-2 vanishes at a stagnation point and that f^(0) = 0, the mass flow rate 
through the wall becomes 


(pv) = - — 4 — — 

w r x v^ 


For the axisymmetric stagnation point, £ may be evaluated as 


/* 2 / du e \ 

-L v.v 41 * v^var), 


4 

X 

T 


where r is replaced by x and u g 



(15) 


(16) 
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The wall -boundary condition on the stream function f is expressed in terms of 
the suction or injection mass flow rate using equations (15) and (16) as 


f 

w 


= f(0) = 


(pv) w 


/du \ ~ 

1/2 

^w^' w\dx /g 



(17) 


Equation (17) may be used to obtain f(0) for given free- stream conditions and wall 
temperature when the injected mass flux (pv) w is specified. 


Method of Solution 

The transformed boundary-layer equations (10) and (11) were solved subject to 
the listed constant conditions. Equation (17) was then used to relate the wall mass 
flux to the boundary-layer solutions. Equations (10) and (11) with the listed boundary 
conditions constitute a two-point boundary value problem. The solution was obtained 
by treating the problem as an initial value problem as follows: values of f^(0) and 

g (0) were assumed; the transformed boundary- layer equations were integrated from 
” 1 

the wall to the boundary -layer edge on a digital computer using a fourth order Runge- 
Kutta scheme with an integration step size of 0. 01; the values f^(0) and g^(0) were 

automatically iterated until the following conditions were satisfied. 


g(00) = l 

g ?? ( 00 ) = 0 

For all calculations presented in this work, the above conditions were satisfied within 
-4 

tolerances of 5 x 10 . A listing of the boundary-layer program is given in the appen- 

dix. 


V“) = 1 


f (” ) = 0 
7777' ' 


*The present problem was programed for digital computation by Mr. Herman 
Hines of the Computation and Analysis Division of NASA Manned Spacecraft Center. 
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Pg 

Correlations of — and l 
P 


Pg 

It has been pointed out that — and l must be functions of r\ or constants in 

order to treat equations (10) and (11) as simultaneous ordinary nonlinear differential 

Pg 

equations. Calculations have been made of the quantities — and ! as a function of g 

(using the equilibrium air properties of ref . 3 and Sutherland’ s viscosity law) for 
higher altitudes and velocities than were considered in reference 1. Velocities to 
30 000 feet per second and altitudes to 220 000 feet were considered to conform more 
closely to flight conditions experienced by manned spacecraft. The calculations shown 
in figures 1(a) and 1(b) indicate that for a given flight velocity a single curve can be 
used to represent each of these functions over a wide altitude range with little error. 

Pg 

It is stated in reference 1 that moderate changes in the distribution of I and — for 

r 

fixed end values resulted in negligible changes in heat transfer. The correlations of 
figures 1(a) and 1(b) were fitted to equations of the form 



0 .2 .4 .6 .8 1.0 0 .2 .4 .6 .8 1.0 

9 9 


(a) = 25 000 ft/sec. (b) V^ = 30 000 ft /sec. 

Figure 1. - Density- viscosity and density ratio correlations. 
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(19) 


^=1- rj(l-g)- y 2 (l-g) 4 

These are identical in form to the equations used in reference 1 when = 0. Values 

of a and y used to fit the correlations are tabulated in figures 1(a) and 1(b). The 
solid lines represent these curve fits. 

RESULTS AND DISCUSSION 


Solutions of the stagnation-point boundary-layer equations have been obtained for 
dissociated-equilibrium air for a range of suction and injection conditions. The 
method employed can also easily provide solutions similar to those of Reshotko and 
Cohen (ref. 4), Solutions of both types are discussed in the following pages. 


Comparison With Results of Fay and Riddell 

A solution has been obtained for comparison with the results of Fay and Riddell 
(ref. 1). This solution corresponds to that listed on the first line of table I of the 
reference. The resulting boundary-layer functions are shown in figure 2 together with 
the input data used, It is evident that the functions vary smoothly and satisfy all the 
required conditions quite well. A comparison of heat transfer and shear parameters 
for this case for L ■ 1. 0 and o« 0. 7 is given in table I. 
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Figure 2, - Equilibrium air boundary- 
layer functions. 


TABLE I. - COMPARISON OF PRESENT 
WORK WITH REFERENCE 1 



Present work 

Reference 1 

g,(°) 

0.2191 

0. 220 

v°) 

.3258 

.326 

Nu 

.3137 

.314 

Vr 




The present program may, therefore, be considered capable of duplicating the equi- 
librium boundary-layer calculations of Fay and Riddell. 

Calculations also were made for the case shown in figure 2 using an integration 
step size of 0. 001 to evaluate the adequacy of the selected 0. 01 step size. No differ- 
ences observable on the scale of figure 2 were found. 
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Comparison With Results of Cohen and Reshotko 


In references 4 and 5, Cohen and Reshotko presented similar solutions for the 
compressible laminar boundary layer under the assumption of linear viscosity- 
temperature relation. These solutions may be computed for the stagnation-point cases 
using the present computer program by taking 


a l = a 2 = 0 

a 3 ~ 

yj =1.0 

to 

II 


Comparison of a cold- wall ^g w = Oj stagnation-point solution with one of the solutions 

tabulated in reference 5 indicated close agreement in the boundary-layer functions f, 
f , f , and g . A few values for this case are listed in the following table to 

7] 7)7] n 

demonstrate this agreement. The stagnation-point boundary-layer solutions in table n 
were computed using the following conditions: 


g(0) = 0, f(0) = 0, a = 1. 0 


TABLE H. - COMPARISON OF PRESENT WORK WITH REFERENCE 5 


Cohen and Reshotko (ref. 5) 

Present work 

V 

f 

f v 

f 

m 

s n 

V 

f 



S V 

0 

0 

0 

0. 5806 

0. 4948 

0 

0 

0 

0. 5812 

0. 4943 

1.0211 

0.2913 

0.550 

.4576 

.4479 

1.02 

0. 2913 

0. 5500 

.4577 

. 4471 

2. 5059 

1. 4738 

.950 

. 1005 

. 1303 

2.51 

1.478 

.9507 

.0995 

. 1287 

3.9916 

2.9384 

.999 

.0032 

.0049 

3.99 

2.937 

.9990 

.0032 

.0049 


The Cohen- Reshotko type solutions have also been generated for a series of cases 
in which varying rates of mass injection and suction have been employed. These solu- 
tions may be used to show the effect of suction or injection on the heat-transfer rate at 
the stagnation point. Convenient parameters for expressing these results are the ratio 
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of the heating rate with wall mass transfer to the heating rate without mass transfer 

a (P v >w 

and the mass-transfer parameter B' = — . Negative values of B' indicate 

«0 <’e u e N St 

mass injection into the boundary layer. Figure 3 shows a comparison of the present 
calculations and results from reference 4 (for a = 0. 7 and g^ = 0). The close agree- 
ment of the calculations indicates that the present program accurately represents the 
solutions of Cohen and Reshotko with mass injection. 



-B’ 

Figure 3. - Effect of wall mass injection 
on stagnation-point heating rate. 



Figure 4. - Equilibrium real-air f 
profiles. ** 


Equilibrium Real-Air Solutions With 
Wall Mass Transfer 

It has been shown that the present 
program accurately reproduces the equi- 
librium real-air stagnation-point results 
of Fay and Riddell and the solutions of 
Cohen and Reshotko for stagnation-point 
mass injection. In the present work, so- 
lutions have been obtained for equilibrium 
real air including the effects of dissocia- 
tion, variable properties, and mass trans- 
fer for flight velocities of 19 700, 25 000, 
and 30 000 feet per second. The two sets 
of solutions for 25 000 and 30 000 feet per 
second were obtained using the density- 
viscosity ratio functions shown in figure 1. 
The three families of solutions are tabu- 
lated in tables in, IV, and V. 

Distributions of f obtained for a 

n 

range of values of f(0) between -0. 6 and 
+0. 6 for the = 25 000 feet per second 

case are shown in figure 4. It is seen that 
the boundary layer becomes thinner for in- 
creasing suction (f (0) positive) and thickens 
for increasing injection (f(0) negative) as 
expected. 

The effect of mass transfer on the 
stagnation-point heat transfer rate for dis- 
sociated air is presented in figure 5. It is 
evident that the calculations for the three 
flight velocities show negligible differences 
when presented in terms of the mass trans- 
fer parameter B'. The independence of 
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2.0 



- 1.6 - 1.2 -.8 -.4 0 .4 .8 1.2 1.6 2.0 2.4 


-B' 

Figure 5. - Effect of wall mass transfer 
on stagnation -point heating rate for 
equilibrium real air. 


flight velocity was indicated in reference 6. 
The correlation curve 


= 1 + 0. 72B' + 0. 13(B’) 2 (20) 

q 0 


presented in reference 6 is seen to be an 
excellent representation of the calculations 
in figure 5. Care must be taken, how- 
ever, in applying this relation for values 
of B’ greater than those for which 
boundary -layer solutions are available. 


For comparison purposes, the Cohen and Reshotko solutions of figure 3 are re- 
produced on figure 5. It is evident that a significant difference exists between the 
present solutions and the solutions of Marvin and Pope on the one hand and the Cohen 
and Reshotko solutions on the other. This is attributable to differences between these 
two classes of solutions in property variation across the boundary layer. It appears 
that the major effect is due to the variation of the density -viscosity ratio £ . Exami- 


nation of equations (10) and (11) reveals that terms involving £ 


and 


d£ 

dg 


appear in 


both equations. The term is important in determining the character of the solu- 
tions particularly in the no-injection case. Equations (10) and (11), respectively, may 
be rewritten for the case d = 0 as 


W7 


1 

£ 



_S_ f 2 +f U 

p r\ rjri dg & rj 


( 21 ) 


and 


s m 


a 

£ 


fg + lg 2 ^ 

\ & T] (7*7 J dg J 


( 22 ) 


d£ 

It is evident from figure 1 that -3- is large and negative in the vicinity of the wall for 
the real air case. On the other hand, ^ vanishes in the solutions of references 4 and 
5. The importance of this term is particularly evident for the no mass injection case 


dl 

at the wall in which g (0) vanishes with -3-. 

°7]7] dg 


The difference in character of the 


no-injection solutions is readily apparent in comparing the monotone-decreasing 
values of f and g of reference 5 (see table n) with the corresponding curves of 


m 


v 
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figure 2 which exhibit maximums in the boundary layer. For cases with mass injec- 
tion, the constant i solutions exhibit maximums in the boundary layer, but the 
maximums are much flatter than for the varying l solutions. Further comparison 
of the table and figure 2 shows that the boundary-layer thickness in terms of 7] is 
much larger for the constant i solutions than for the variable l solutions. Thus, 
in view of these fundamental differences in the variable property solutions of this 
report and in the constant l solutions of references 4 and 5, the difference between 
the blowing curves of figure 5 is not surprising. 


Injectants Other Than Air 

The present solutions of the boundary- layer equations have been obtained for air 
and apply strictly to the injection of air into an air boundary layer. It is likely 
that these results will also be quite accurate for injectants whose diffusion and thermal 
properties are similar to air. Theoretical and experimental studies have shown that 
greater reductions in heat transfer at a given value of the mass-transfer parameter 
are obtained when the injected gas has a low molecular weight or a high-heat capacity. 
Stewart (ref. 7) examined the experimental and theoretical information available on the 
injection of foreign species into a flat-plate laminar-air boundary layer in the absence 
of dissociation effects. Foreign species considered were gaseous water, helium, hy- 
drogen, and carbon dioxide. It was found that the results for the above species could 

be approximately correlated in terms of the parameters. and B' ^ c p*y ^ where 

c p * is the ratio of the specific heat of the injected species to that of air. The specific 

heat- ratio values of c *, as obtained from reference 7, are given in the following 

P 

tabulation. 


Air 1.0 

Helium 5.26 

Hydrogen 14. 3 

Water vapor 1.975 


Carbon dioxide . 0.854 

0 4 

If it is assumed that a similar correlating parameter B* ’ applies to the 

stagnation-point boundary layer, the preceding table may be used to estimate the 
reduction in heat transfer in a real-air boundary layer as a consequence of the injec- 
tion of an arbitrary foreign gas. This is obviously accomplished by simply interpre- 
ting the abscissa of figure 5 as B’ 4 - This approach is probably most accurate 

at low and moderate injection rates and low wall temperatures where real-gas effects 
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involving the injected species may be neglected. More recent studies (for instance, 
refs. 8 to 10) suggest that simple correlations, such as that given in reference 7, are 
inadequate to describe fully the effects of foreign-gas injection. Therefore, caution 
should be exercised in the use of the suggested approximation. 


CONCLUDING REMARKS 


The equilibrium -air stagnation-point boundary layer with suction and injection 
has been examined with the use of a computer program. The program inputs allow 
calculations for a range of assumptions regarding the property variations through the 
boundary layer. It has been shown, for instance, that the solutions of Reshotko and 
Cohen and of Fay and Riddell may be generated. Solutions may also be obtained for 
those gases other than air whose properties may be expressed in terms of the selected 
property-fit parameters. 

Calculations of the real-air boundary layer with mass injection have been found 
to be in good agreement with the correlation equation of Marvin and Pope. The real- 
air mass -injection curve of the present report is found to exceed the corresponding 
curve obtained from the low -temperature solutions of Cohen and Reshotko. The differ- 
ences are attributable to the different forms of variation of the density -viscosity ratio 
across the boundary layers. 

The effects of foreign-gas injection on stagnation-point heating may be estimated 
from the present results by the use of a modified mass transfer parameter involving 
the specific heat of the foreign gas. 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, March 26, 1968 
124-07-01-01-72 
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TABLE in. - REAL-AIR STAGNATION- POINT SOLUTIONS, V = 19 700 FT/SEC - Continued 
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TABLE JH. - REAL-AIR STAGNATION-POINT SOLUTIONS, V M = 19 700 FT/SEC - Continued 
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TABLE m.- REAL-AIR STAGNATION- POINT SOLUTIONS, = 19 700 FT/SEC - Concluded 
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TABLE IV. - REAL-AIR STAGNATION-POINT SOLUTIONS, = 25 000 FT/SEC 
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TABLE IV. - REAL-AIR STAGNATION- POINT SOLUTIONS, V^ = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION -POINT SOLUTIONS, - 25 000 FT/SEC - Continued 
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TABLE IV. - REAL- AIR STAGNATION -POINT SOLUTIONS, V = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL- AIR STAGNATION-POINT SOLUTIONS, V w = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION-POINT SOLUTIONS, V x = 25 000 FT/SEC - Continue^ 
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TABLE IV. - REAL-AIR STAGNATION- POINT SOLUTIONS, * 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION- POINT SOLUTIONS, V M = 25 000 FT/SEG - Continued 
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TABLE IV. - REAL-AIR STAGNATION- PCttNT SOLUTIONS, V = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION- POINT SOLUTIONS, = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION- POINT SOLUTIONS, = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION-POINT SOLUTIONS, = 25 000 FT/SEC - Continued 
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TABLE IV. - REAL-AIR STAGNATION- POINT SOLUTIONS, = 25 000 FT/SEC - Continued 
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TABLE V.- REAL~AI$ STAGNATION-POINT SOLUTIONS, * 30 000 FT/SEC 
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TABLE V. - REAL- AIR STAGNATION-POINT SOLUTIONS, V*. = 30 000 FT/SEC - Continued 
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TABLE V. - REAL-AIR STAGNATION- POINT SOLUTIONS, = 30 000 FT/SEC - Continued 


UI UJ LU UJ UJ 
O' ^ 

CM O' .h 

O'tOCO't 
*0 O' ^ 
mm® cm n- 
O' 'O *h N- TO 

ino if* O' 

O’ O' co n- n- 


O min in 

CM N- <30 4> f\| 
'OONom 
o o >o o 

co m imn in 

^ (Min-^o 

O' 4 O' in pH 

vo 43 in in in 


UJ UJ LU UJ IXJ 

oinomco 
rg 4 <n cm in 
(ONMN4 
43 43 o m o 
to o in in o 
niN'f't h 
r- m o f- 4 
4 4 4 m m 


uj uj ui uj ui 
'O O' Ol 'O CD 
4) 00 CO O' O 
4- N- 4- O fO 

-o in cn in n 

N- >0 in (\J 43 

H E N CO O 

cm O' r- in 4 
m og cm cm cM 


UJ UJ UJ UJ UJ 

in in oo O' oo 

co m to co m 

4 0'co»hn. 
n- to in rn 4> 

43 r-» O' o cm 

4- O 43 in 4 

CM rH O' co n- 

PO N r4 H rt 


UJ UJ UJ UJ UJ 

cm m w* n- 4 -t 

IM ffl O H N 

CN n- o> on 4 

N® HON 

in OD H (\J o 
4 in cO <— i in 

4 in 4 4 ® 


4 H CN 4 CM 
h 4r m cm »h 
n ® o m 4 
OOiENH 
O' 4 O N- 4" 

NNNHh 
H H H r< 
• * * ♦ • 

OdOOO 


UL* UJ IX* UJ UJ 
4 fM 43 CO O' 
C-» o 43 in CO 

r-n oo 4 co rg cm 

43 o oo m H 

o o m m o 
NNN® o 
i m in n- o 


UJ UJ ui UJ UJ 

cm 4 in oo O' 
4 in o 4 o 
4 O N 4 4 
n- 4 cn o to 

4® h4 h 

cm in o in cm 
CM 4 N- O' eg 
CM CM CM CM CO 


UJ UJ UJ UJ UJ 
O CM 00 4) 00 
4 m co m rg 
4 O' 4 O' CM 
•H 4 O O' .O' 
00 43 4) in 4 
O' CO co O' —• 
4 N* O TO N- 
m cn 4 4 4 


UJ Uj UJ UJ UJ 
HNHNin 
pn r- CO CM cn 
CM 4 O' CM in 
N 4 h o ® 
o o o in co 
4 r- o cm to 
OWNOffi 
in m m 43 43 


UJ UJ UJ UJ II i 

—i co m a in 
h rn in 4 m 
O cn 43 N* o 
in n- 4 o oo 

4 m N- 43 O' 
m o m cm in 
43 O' ^ TO 4 
43 4 N N.N 


UJ UJ UJ UJ UJ 
43 O O' 4 in 
4 © rn in O' 
in o cn o> 4 

43 in 4 ,#h 00 
O' N- N* 43 cn 
CM eg 4 cc 4 
in in 4 cm o 
N- N- !h- N- N* 


UJ UJ UJ UJ UJ 
.-I 09 O' CM CM 

O' cm O' cn cn 
minnino 

rH in O CM 00 

4 4 43 CM 00 
CM cn 00 O' 43 
n- cn co cn co 
4> 43 in in 4 


o c> o o o 

lilt 1 

UJ UJ UJ UJ UJ 
o n N O' cn 
o m O' m 43 

OIMENN 
o O N- co 4 
in o in 4 o 
4 co — < in o 
43 CM O to to 
H CM cn m 4 


o o o o o 

fill 

UJ LLI U ‘ UJ UJ 

N- in in cm »— i 
cn ,-h 4 o o 
4 4 cn m m 
to r- m 4- 4 
43 43 43 0 N- 
tn fM — t en rg 
co o cn n- -* 

IA N C O' H 


UJ UJ UJ UJ UJ 

O' O «"-* CM 4 

co n- in in 4 

'Onontn 
m 'J- CM CM O' 

in n- 4 4> m 

O' N- N- 00 ' 

CM 4 4) CO »H 

•— i •— I H h M 


LU U ‘ UJ LLI UJ 

m 4 to in o 

CM CM N- O N* 
cm O' go m co 
co o co o H 
N- CO 4 CO N* 
in w O' co o» 
m 43 00 H 4 
cm cm cm m m 


UJ UJ UJ UJ UI 

O' o in iH 4 

CM N- 4 •— I N 
H O' O' O' ® 
'OHifflrlO 
O HfO ® 
CM 43 H N* 4 
co •— i in co nj 
m 4 4 4 in 


LU UJ U3 UJ UJ 

m O' in 43 in 
O' h co O' cn 
o in in m cm 

cc in CM 00 
CO 4 O' CC 41 
<-* O' 43 pn O' 
43 O' m N- o 
in in 43 n- 


UJ LU UJ UJ UJ 

rn O O' 4 4 

CM e CM 43 CM 

m n- o «-i cn 

m o ® O' in 

to m in rg 

4 4> N- in •— i 

4 |N o cn 43 

r- n» cc co co 


4 N- m 4 o 

in 0 to in o' 
4 'C Cj N* to 
H ® 4 m N- 
0* -i N 4) O' 
CMOtoCMT- 
N* N- 4) 4> in 


ociooo 
o o o o o 


o o o o 

I, i I I 

ui LU UJ UJ 
P rH 43 O 

O- © CL. (M 

©to 4 ■*-• 
© O o 

H ^ O' N 
4 m N- N- 
o c\j in © 
*-i cm m to 


LL. UJ UI UI UJ 

o cc in o m 
Q © 43 N- O 
C 43 m 4 > o 

.o 4* cn a cm 
o n- a 4 m 
O O' cn f- it. 
O O' O' O' o 


0 c o o o 

1 I; I I I 


LU UJ U LU 

O O' O O 

cj cr c- o 
Cl a- c o 
cJ cr c c 
dcoo 
cj cr c c 
cJ O' in o 


UJ UJ IU U* UJ 
rq 4 4 O' c6 
O] 4 to 4 on 
4 h m a- h 
to m co cc co 
cd in o' rn O' 
rg m 4 43 co 
IN in o cm 4 

4j d' r—J »*”< . 


tu uj UI LU UJ 
in 4 43 m oc 
o n 4> O' eg 
4 4 CM O' N 
CCj N- CM 4 N~ 

cn in co O «-< 

CM ce m to -I 
oj CC CC N N 

*—j ,~l !— < r-j rH 


P C. C O o 

C C pi o c 

UJ UJ UJ Ui UJ 

0 C > CN o O 

01 C. O' O c. 
O. o O' o o 
C O O' o o 
O O O' o o 

CJ C O’ O CJ 
IP O 4 C IT 
c\ m m 4 4 

• • • • • 

O C_> o Q O 


, UJ Ui U. UI Ui ! 

I o 43 4) 43 N f 

cn .43 in oo oo j 
. 4 43 N- m 4) t 

j 43 cn to O' ! 

i cn m a, cn 43 i 

! cn O' 43 in <1 : 

i N- CN CM in CO I 

I H h (N M N 


; lu uj u uj uj j 

! 4) 4 N- 4 CO I 

i cm m cm 43 in | 

i in 4 in o o j 

] CM 43 CN 4 i 


UI U UJ UJ uj 
43 O O 4 
O' 4 cn N- P 
CD O' ♦“* 43 4 

O' in O' in co 
o w-4 to r- in 
to CM 43 CM O'. 

^ ui re cm in 

m m m 4 4; 

* • • • i 

o o o o o 


o o o o c? 

Ill I 
UJ UJ ULJ UJ LU 
CM CM 43 to Ol 
43 O 43 to 4 
m N- 4 4) N-; 
•H N- 4) O' 4l 
g 4 N cr od 
CM 43 O' to; 
O to N- 43 4 
*-♦ to 43 4 nj 


43 N- 4 41 to 

in in —< to cm 

43 to N- Oj CO 

.H CM to H to 

h- to m H in 

o cn n- i-i 4 

4 to to 4j 4) 

• • • 1 4 


m w~i c- 
ooocjc 
I I I I, 

UJ Ui UJ li| u> 

t m to m n- 
N- CM O in O 
m cn 4 ) oi cm 
in O O toj CT. 

4 4 4 03 N 
to to to -H m 
4/ 4 CM CM] -4 
O m to. cd —• 


3 O O Cj c 


LU UJ ui UJ UJ 

43 o to to «-< 

H IN (J H O' 

43 cm m m 43 

or-moN- 

to to to O p-4 
4} CM to CN O' 
TO CM to CO pH 
43 N- IN N- CO 

• • 4 * * 

o o a o o 


rg tr Nt ^-t m 
TO O' CM N- O' 
o o CM pH to 
4 4 Ifj C H 
»h m m o cm 

N (\i O CD to 
4 CC <- to CN 
pH ^ cv CM rg 


Uj LU U UJ UI 

a o c o o 
a o c o a 
a o c o a 
o o c o a 
o o c o o 
to o if o in 
cm m ce 4 4 


1 c c c do 


68 




TABLE V. - HEAL- AIR STAGNATION- POINT SOLUTIONS, V w = 30 000 FT/SEC - Continued 


I 

■ y 


eq 

9 


Si 



O 

I 11 




Q 



o o o o 


LU UJ LU LU UJ 

irt in — < 4 «-< 

O <0 ^ O' m 
N >o O' in 
4 w >o cm in 
in w 4 O' -O 
•-4 © N- in 4 
•-<0000 


O O O O O 


O O O O O 

O o o o o 


Hi UJ UJ LU LU 

p in © tn cm 

N N © © 4 
(O ® O >0 N 
O' 4 00 O' c\j 
HMJ'OO 
(*) O'N C"4 - 
d n cm n p 
4 p p cm cm 

G> O O O O 


o o o o o 

O O O O O 

LL) LU ID LU LU 

h* in in o> 4 

n « H CO 4 

n 4> oiMm 
m 43 cm p 
cm in csj 4 cm 
< r 4 * cm n o 

® O N (nin 
CO O' .0* O' O' 

o o o oo 


§ 


o o o o 

o o o o 


UJ U! UJ UJ LU 

© O' -4 n oo 

®H<U)0' 


3 


■O cm O' CO 
cx> <r o cm 
© p © n 
'O o o 

h 'f H 

P\1 • — * »— I ri 


Q O O O O 


Ui LU UJ UJ 
CO WOO o 
4 CM CM 4 P 
*4 O' HO O' 
tv 4 in w 
n tv CO' n 

If IV h 4 H 

4 in o mo 
O' O' O' O' 

• *» » 

©OOP 


q © o < 

ui UJ U l 

4 'O p c 
d iv pi , 

cm p to > 

dN C0' 
p © © L 
ui cm p ( 
q <4 < 

uj-OiOt 

q © 4 < 


^ h 4 H d 

c o o o q 

U UJ U| LU ui 

C O Q O 0 

c © q o © 

c o d o q 

c o o o d 

c © a o a 

I v in 



o o o o o 

UJ UJ LU LU LU 
4 P M3 O CM 
CV 00 4 S tft 
CM TO 4* O CM 
CM CO CM O O' 
ifl N H 4 H 
rfllVNHH 

o o o o o 


o o o o o 


0-0 o o »-* 
o o o o o 

i > 

UJ Hi UJ Ui UJ 
OHOMCI 
4 CON INI H 
CM O' O O oo 

® on •t in 
co in 4 r-< o' 

nNOOH 

O' in n o oo 

►— ( rH »— < r— t N 
. • * • • < 

o o o o c 


UJ LU UJ LU U. 

m hn Ntn 

HNfflHO 
O' H O N O 
O' in P T- —* 

oc «o N p oo 

O0"0fVi0 
>© 'O N- TO 00 
© O' O' tf O' 

• 4 4 4 • 

o o o o o 


00 oo o 

1 I I I I 

LU UJ UJ UJ LL 

-4 4 - in p © 

O 00 N 00 © 
© 4- 43 P 

© 4- in >0 sC 

03 © cm o in 

in o n © o 

tn in or in 43 

TO 4) 4 P (N 

• * * * * 

O O © O O 


OO o o o 
o o o o o 

UJ UJ U I UJ LU 

P 4} r-t © 4 

4- «-t 1(/ N 43 

N- 43 in N m 

O in co co tv 

h oo lO Ntn 
4) © CM 4* 43 

OC 00 © © © 

© © © © © 

# • * * • 

o o o a o 


o © o o -< 

o o o o o 

UJ UJ UJ LU LL 

<G h 4 (V N 

O if r CD if 

in p if' 4 * cd 

4" 

4) 
r- 

© in a in o 

ffi «d O' O' H 

• t • « • 

© q o a o 


4n o c 
o 4> p o 
r- 43 4} o 


HrlH 

o d o a o 

LU LU UJ UJ UJ 

o o o o o 
o p o Cl o 
o op a o 
o o © a o 

o o © O o 
o m o in o 

OOHHN 

CM CM CM CM fM 

• * • * « 

o q o a o 



o o o o q 


UJ UJ LU UJ UJ 
min-tHO 

-i oo © p cm 
oo co co in n 
r~- co in cm cm 

TO 43 4 P CM 

o o o o o 

o o o o o 

hhhhh 
p O O O P 



LU LU UJ UJ LU 

n- co © tn -< 

cm 4 m m p 

43 in © 4 N* 

in O 43 4" CM 

f~< O P P 

o o o p d 

p o o o o 

HHHHH 

P O O p P 


t-4 r-t O O 
O O O O P 

UJ UJ UJ UJ LU 
O'NO'OH 

in N- in © 43 
in N (VI O 43 
—»0 0©N- 
O O O © © 
O O O © © 

o o o © © 

**H r-A © © 
• • * * * 

O O O O O 


O O P o o 
11)1} 
UJ UJ UJ LL UJ 

® h 4 in a 
H N ® P 4) 
N N in N Nr 
© N- in © © 

CM 43 ■*-* TO p 

nn tn P in 
o tn 4 in oo 
4) 4 P N H 

O O P o o 


-H CM CM CM CM 

0 P © O O 

1 J I I J 

UJ UJ LU LL UJ 
43PN-pN> 
— < 4 43 P © 

03 CM c in N 

4 © a P 
43 m 4 N> 4} 

P P N- CM 4> 

P in 43 sQ rH 

>-< © 43 4* P 

• • • • • 

OO OOP 


CM CM P P P 

0 O O O O 

1 I l i i 

LU UJ LL LU UJ 

P © P © 4 
© 4 O P -i 
©4^04 
O 43 4 O 4 
© i0 P iO 43 
P CM © O CO 
4 P ^ 4) 
IN) H O' >0 4 

• » I « • 

0.0 POO 


o o P o o 

oo o oq 

LU LU UJ LU LU 

cm 4 © cm in 
o p oo © p 
P © © CM P 
43 © 03 N- 43 
CM 03 03 P 4 
O nj.4 43 N* 
© © © © © 
© © © © © 
• • • • • 
O O o O o 


o o o o o 
o o © a © 

UJ UJ UJ UJ UJ 
43 in © P CO 
43 4 *“t p tn 
4 r-t (N •— I 4 
m CM CM O © 
CM P CM in 43 
P P © © © 
© © © © © 
© © © © © 
• • • 4 • 

o o o q o 


o o o .«-« -i 

o o o o o 

UJ LU LU UJ UJ 
min O' him 
43 43 © O 43 

m m h h (v 

nh n p o 
» © © o o 
© © © o o 
© © © o o 

© © © rH *— I 

• • 4 4 4 

© o o o o 


h h CM CM rsi 

0 o o o a 

1 i i i t 

LU UJ LU LU UJ 
m CM ©ON 
00 N cm.n 4 
4 h 4 Q iR 
4 P 4 ^ 4 
CM w 4 t-H 4 
43 -h © r-t N 

ra m a <M r| 

W ~ © 43 4 

d C' ood 


CM CM CM P p 

o o o o o 
1)111 
LL UI UJ LU Ui 
N in cn CM P 
4 P CM Nj N 
P IT, O P P 
4 P -< 4 P 
n in © cm N 

in p cm H P 

N N w Of O 
CM H H 4J 4 

ooodo 


m44 4 4 
o a o o o 
i i i i i 
LU Ut UJ UJ LU 

N © CM 4 •-* 

O O P O N 

H H fr N P 

43 in N 0 4 

CM CM r- 4 © 

43 4 in p a 

t-H 43 a 4) 43 

CM © CM CM 4 

4 4 4 a 4 

o o q o o 

i i 


LU UJ LU LU U 

4 O P 43 If 

P O N CM a 

p O © P ^ 

P CM P 4 C 

4 {\i N *4 

N C W O' O' 

© LT U' © O' 

© © © O' 

4 4 4 * 

o o o o q 


LL UJ UJ lU I 


CSJ . 


i in i 


© © P P 
IT N © 4 
md in q. 
in 43 r- p 
© © © d 
© © © d 
© © © oi 

o p o q < 


UI UJ U| Ui UJ 

nioot on 

m ® in m h 
43 P <1 43 in 
4 in <n m in 
p co P TO P 
© © © © © 
© © © © O' 
cr (j> U'. 

> • 

o o q o o 


o o o o c 

UJ ui UI UI u 

N'CNNC 

o CM P 4 ) C 

p CM 4 © 4 

© 00 N 43 4 

p in in in ir 

m o in c if 

Q H H CM (V 

p d *— 1 •“ 

4 4 4 *4 

d o q o c 


liod 


docji 

ui uj ul iu u 
o o o o c 
q o o c c 
d O C! o c 
a o q o c 
o.o a o c 
in o in O' tr 
cm m p 4 
cm cm P cm cv 

J • * 4 

q o q o < 


O O p c 

LL Uj uj ui l 

n tn N a co 
in N p cs p 

pf 4-4.0..© < 
43 43 4 if l 
tn in tn if i 
oinoir _ 

p P 4 4 tn 

HHHr t 

• 4 

o o o d o 


; hJh, 

; o q O I 

LU LLj UI I 
O Ci O 

o a o © 
o o o © 
o a o © 

OOP© 

o in o 

m Ifi 4? 

. CM CM CM 


in oc N p © 


N 43 43 

in it. irj 
in in in 
in o in 
m q 


in 4 
in in 
in in 
o p 
N N 


oo # qo a 


O C' c 
o o c 
o o c 


© o 
© o 
© o 


o c q © o 

- -- © C 4 ji 


69 




TABLE V. - REAL-AIR STAGNATION -POINT SOLUTIONS, V = 30 000 FT/SEC - Continued 
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2200000 OF 01 -0.394735566 00 0,30^69417^ OC 0.29237233F 00 0. 170 5692 6E 00 0.2200 33 82 E 00 0.39974552E 01 
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3500000 OP 01 0.36210245F 00 0.B9868421E 00 0.363749426 00 0.804700916 00 0.55013774E 00 0.121156766 01 
35 5000006 01. 0.40 747 4..? 7 6 00 0 .9 15J 7632 E. 0.0..„. . 0.3238 Q220E 00 0.8312 54 40 E 00 0.51127001E 00 0.11776799E 01 
3600000 0E 01 0 .45 3656 56 E 00 0.931060896 00 0.28364556E 00 0.85576472E 00 0.46861970E 00 0.11480449E 01 
365 00 00 QE 0.1 0 , 50 0 54 76 6E 00 0 • 9.4 425 648E. 0.0 . . . . _G... 2.4 4425.18.6 .0(1. .Q.87807506E 00 0.42345698E 00 0.112234Q5E 01 
37 000 00 06 01 0 .5 480 502 2E 00 0.95553596E 00 0.2071 5529E 00 0.89809588E 00 0.37726539E 00 0.11002361E 01 




TABLE V. - REAL-AIR STAGNATION- POINT SOLUTIONS, = 30 000 FT/SEC - Concluded 
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4349 999 9 E 01 0 • 16900512E 01 0.10000417E 01 0. 8535 664 8E-Q4 0.99998532E 00 0 . 57 559478 E- 03 0.100001 3 IE 01 

43|7C00OOE 01 0. 17 100 520 E 01 0.1 00004 33E 01 0. 7273 78 5 6E- 04 0.99999587E 00 0 .48 352774E-03 0.1 000 003 75 01 




APPENDIX 


COMPUTER LISTING OF THE BOUNDARY-LAYER PROGRAM 


SIBJOB COBLP GO 
S1UFTC COBLP 

DIMENSION F<3) >Y<3) >TP(2)fFK(3) »AY(3) »AYP<3) »SUM1{2) »SUM2(3) »G0G(4 
.000) *AF0125) »AGP0(25) »AFDP0(25) fFNCT(lOOO) rFP(lOOO) rFDP(lOOO) »G(10 
.00) »GP( 1000) » ETA (1000) »RB( 1000) 

. »XXU2)»YY(12) 

READ(5»200) IPRINT 
200 FORMAT (13) 

INDX=0 

W'REAU(5»100> NOFV » IPSION 

10U FORMAT 1213) 

ISDT=U 

IF(NOFV.EQ.O) GO TO RO 
ICTK=0 

READ (5> 101) NP»SIG» ALPH1 » ALPH2f BETAl » BET A2» GAMA 1 »GAMA2 »FPO»FDPO 

. . GPO r EF» EG » DELF » DELG ► ENN* TW * PS * RHOS * HS* RFT * PINF * R * LWN * ALPH3 

101 FORMAT <I4/(7F10.0) ) 

KK=U 

IF ( IPSION.NE. 0 ) GO TO 500 
READ(5>102) GO 
GO TO 501 

500 READ (3* 102) (XX ( I ) * 1=1 * IPSION) t ( YY < I ) * 1=1 * IPSION) 

60=XX(1) 

301 CONI INUE 

INDX= INDX+1 

GSAVE=GPO 

FSAVE=FDPQ 

READ(3*102) (AFO(I) *I=l*NOFV) 

102 FORMAT (7F10.0) . 

111 ICTK=ICTR+1 

ISP 1-0 

FOSAFOIICTR) 

IFtFO.NE.O.) 60 TO 103 
GPO=GSAVE 

FDPO=FSA\/E 

GO TO 10R 

103 GPO=AGPO(IFG) 

FDPOSAFDPOIIFG) 

GO 10 10R 

4TT READ (5*1) NP * SI'G * ALPH1 * ALPH2 * BETAl ♦ BETAS - * GAMA 1 * SamASTB 7FP0*FDP0* 

. GPO * EF * EG * DELF * DELG * ENN * TW * PS » RHOS *HS * RFT , P I NF * R * LWN ♦ ALPH3 
KK=0 

1 FORMAT(I4/(7F10.0) > 

READ (3*6) FO 
6 FORMAT CF10.0) 

IF (IPSION lNE.'”0) 53 TD“b02 

READ (5# 102) GO 
GO TO 503"' 

302 READ(5»102) ( XX( I ) * 1=1 * IPSION) * ( YY( J) *J=1* IPSION) 

G0=XX(1) 

303 CONTINUE 

INDX=INDX+1 — - — 

10R WRITE(6*90) ALPHl * ALPH2 * BET A1 * BET A2 * GAMA 1 * GAM A2 ► B *GPO*FO*FPO*FDPO 
. , DELF ► DELG * EF * EG * ENN * TW * PS * RHOS * HS * RFT * PIN? * SlG * R , LWN * ALPH3 
^0 FORMAT ( 1H1 * 
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. 2X.18H ALPHA-SUB-1 »18H ALPHA-SUB-2 t 18H BETA-S 

• UB-1 * 18H BETA-SUB-2 »13H GAMMA-SUR-1 »18H GAMMA- 

• SUB-2 /bE18.8//2X»18H B ' >18H 61 _ ~»'l 

.OH F * 18H FI . 18H F2 /5 

.E18.d//2X.18H DELTA F2 »18H DELTA G1 »18H TOLERANC 

.£ ON FI . 18H TOLERANCE ON G *18H TOLERANCE ON Gl /5E18 . 8//2X r 1 

.HH WALL TEMPERATURE * 18H P-SUB-S *18H RHO-SUB-S »1 

• OH H-SUB-S * 18H RADIUS /bE18.b//2X» 18H P-SUB-IN 

.FINXTY » 18H PRANDTL NUMBER 18H GAS CONSTANT .18H LEWIS N 

.UMBER • 18H ALHPA— SUB-3 /5E18.8) 

IF ( 1PS10N.NE. 0) GO 10 b05 
«.RIIE(b.b0b) GO 

BOfa FORMAT 1//2X.1HH G /E18.8) 

GO 10 b07 

bUb »vRI TE(6»b08) IXX( I) # YY( I) » 1 = 1. IPSION)" 

bU8 FORMAT (//2X.18H G » 1811 0 /2E18.P) 

bU t CONTINUE 
1FG=U 

wRl lElb.llb) 
lib FORMA I (2X»//) 

CALL CLOCK TTYME) ' " ~ : . .. T „ , 

Ml IH 6 .2000) IT ME 

2UUU FORMAT (12H CLOCK READ F7.2.39H WHEN INPUT PARAMETERS WERE INITIALI 
./ED//> 
lTI>iE-0 
GO 10 2004 
b 0 CALL CLOCK ( T V ME) 

I TIME -IMME+I 
y.RIIt<b.20U2) TY.1E » 1 1 1 ME 

2UU2 FORMAT ( 12H CLOCK REAU F7.2.28H AT END OF ITERATION NUMBER 14//) 
2UU4 1ST=0 

1FG=1FG+1 

AFUPOt IFG)=EDPO 

AGPO ( IFG)=GPO 

ISDI=lbDT+l 

1NDU=J 

ICT=0 

It 10 1 —0 

II = 0 
IFP=U 
IG=U 

12 X=0. 

NEQ-3 

NS0=2 

Hr. 01 

YP ( 2 ) r.GPO 

Y (2) =30 

YP(1)=FDP0 ■ 

Y ( .1 ) =F PO 

Y ( 3 ) =F 0 

JO 0 1=1 #NP 
IF(l.NE.l) GO TO 2 
1ND=4 
GO 10 4 
2 INU=1 

b CALL RUNKUT (Ni_Q* NSO.H. X* F » Y » YP i FK » AY » AYP» 5U-*l »SU ,2. IN".) 

4 IF ( T ( 2 ) .GE.O. ) GO TO 175 
..HI lEtb.l/b) 

1Tb FORMAT l 39H 6 IS NEGATIVE. CASE y. I LL NOf CONVFRGi . ) 

IF (NOFV . EU . U ) GO TO BO 
I F (1 C TR.EQ. NOFV > 60 I 0 99 
GO 10 111 

1Tb SKG=SNHT(Y<2) ) 

GINVSI •/ ( Y (2 ) ) 
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GINV2— GINV*GINV 

FLOG= ( ALPH1*SRG-ALPH2 ) +GINV+ALPH3 
GP0G2=YP<2)*GINV2 

P1_0G=- ( . 5*ALPH1*5RG-AI_PH2 ) *GP0G2 
EB2 = EXP (-BETA2/Y 12 ) ) 

If ( IPS ION .N S. 0 ) GO TO 510 

D0G-BETA1*EB2 

DPOG= IBETA1*BETA2*EB2 ) *GPOG2 
GO TO 512' 

51U CONI INUE 

DO ^00 M=2»IPSI0N 

N-M ______ _ 

IF(Y(2).LT. XX (N) ) GO TO 90 _ 9 
yuo CONTINUE __ _ 

‘909 DPOG= ( Y YTn ) - YY (N- 1 ))/ ( X X ( N) - X X ( N- 1 ) ) 

Q06 = UP0G* (YI2)-XX(N) ) +YY ( N) 

512 CONI INUE 

C0G=-1./IFL0G*U.+00G) ) 

FF~=Y<3> 

P=YU> 

SQP = Y (1)*Y (X J ' 

PP=YPU) 

GM1=1.-Y12) 

GM15=(GM1)**3 

RHOB= X . - ( GAM A 1 +G AMA2*GM1 3 ) *GM 1 

F ( 1 )— l — l • /FLOG) * I (PLOG+FF ) *PP+8 *(RH08-SQP) ) 

F ( 2) =COG* CPLOG* Fl.+DOG) +DPOG*FLOG+sTG*FF ) *YP'<2) 
F(3)=YU) 

IF(II.EQ.X) GO TO 5l 

IF ( Y ( 1 ) .GE..99) GO TO 52 

jpp = iFp+i ' ""~ r ' 

52 IFCY12) .GE..99) GO TO 53 
IG=XG+1 

53 IST=IST+1 

GOGTrSTT=YTFJ 

51 IND=IND+X 

IFUNU.LE.H) GO TO 5 

IF ( I NUT «NE«1J GO TO 22 ___ 

’ NU=2 

GO 10 14 

~ 22 IFUNDT.NE.2) GO TO 16 
MU=3 

GO TO 14 . 

lb NU-1 

INDU=INDU+r ‘ 

R0T2=SQRT(2. ) 

GWALL=1.-G0 

RENU= ( R0T2/GWALL) *GPO 

RH0bU=l./KH0B 

ETAIIJ=X • 

FNC r ( l ) = Y ( 3 ) '■ 

FPU)=Y(1) 

FQpn f=Yp m 

G(I>=Y(2) 

GP ) -YP T2) 

RBU)=RHOBB 
GO TO 3 

14 IF( AB5UPT2)) .LE.ENN) GO TO 7 
3 IF(ABS(YP(2J ) .GT. ENNJ GO TO 201 
FTE5T=ABS(1.-YT1) ) 

”GTE5T=ABM1.-Y(2T) ’ 

IFO-TEST.GT.EF) GO TO (10.15.17) »NU 
IF (GTEST.GT .EG) GO TO UO .15.17 > .MU' 

ICT=0 



DO 150 JSiVINDU 
KK=J-1 

ICT=ICT+1 ‘ 

IF(J.EQ.l) 60 TO 118 
IF(ICT.EQ.l) 60 TO 121 
IFUCT.EQ.5) 60 TO 120 
G0~T0_TT0 
lie WRITE (6» 109) 

109 FORMAT (2X»//) 

WRITE (6 » 116) 

116 FORMAT (2X» 17H STEP NUMBER iTTR F 7T7H FT 

t 17H F2 t 17H 6 >17H 61 

. » 17H RHO/RHO-S/) 

121 WRITE(6>109) 

60 TO 110 
120 ICT=0 

110 WRITE 16»9) KKfFNCTtJ) *FP<J) »FDP(J) r6(J) »6P(J) »RB(J) 

9 FORMAT (2X» 14 » 10Xf6E18.8) 

150 CONTINUE 

SAVE=FDP(1) 

DO 151 K=2*INDU 

IF (5AVE»6E«FDP (K) ) 60 TO 151 

SAVE=FDP(K> 

151 CONTINUE 

IF(SAVE.LE.l.) GO TO 152 

YT=SAVE+.l 

60 TO 153 

152 YT=1.5 

153 XL=0.0 

XR=ETA ( I ) 

Y8=0.0 

WRITE ( 17 1 180 ) FNCT(l) 

180 FORMAT ( 1H1 i40X > 36H THE VALUE OF F(ETA) AT THE WALL IS E18.8) 

CALL 6RID(123>1023f2 t T»924 > XL>XR»YB>YT) 

CALL PL0T1V ( 1 » 1 rETA »FP» INDU* 1 r 1H ) 

CALL PLOT 1 V ( 1 > 1 >ET A » FDP > I NDU rl f 1 H ) 

CALL PLOT IV 1 1 • 1 »ETA r 6» INDU» l.» 1H ) 

CALL PLOT IV ( 1 » 1 r E TA > 6P> INDU» 1 »1H ) 

CALL OUTPUT 

60 TO 75 

201 IF t I PRINT .EQ. 0) 60 TO 8 

KK=K.K+1 

IF(KK.6T. 1) 60 TO 203 

WRITE (6» 116) 

203 WRITE<6» 9) KK rFNCT ( I ) r FP ( I ) * FDP< I ) » 6 ( I > #6P( I ) rRB( I ) 

8 CONTINUE 

7 60 TO ll0fl5»17)»NU 
10 INDT=1 
11=1 

FREF = Y(1) . 

6REF=Y(2) 

FDPS=FDPO 

FDPO=FDPO+DELF 

60 TO 12 

15 INDT =2 

FNEW=Y ( 1 ) 

6NEW=Y(2) 

0ELFN=FNEW“FREF 

DEL6N=6NEW-6REF 

A1=0ELFN/DELF 

A2=DEL6N/DELF 

C1=1.-FREF 

C2=1.-6R£F 

6PS=6P0 
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GPO=GPO+DELG 
FDPO=FDPS 
e0 1Q 12 
17 FNEW=Y(1) 

GNEW=Y<2> 

DELFN=FNEW-FREF 

delgn=gnew-gref 

B1=U£LFN/DELG 

B2=DELGN/DELG 

CELT AF=<C1*B2-C2*B1 ) / ( A1*B2-A2*B1 ) 

DELTAG= ( A1*C2-A2*C1 ) / ( A1*B2-A2*B1 ) 

IND T=5 

FDPNU=F0P5+DELTAF — - — ~ 

GPNUSGPS+DELTAG 

FDPOSFDPNU 

GPO=GPNU 

GO TO 50 
7b GBAR=32.2 

RHOWSPS/ (R*TW ) 

SRTWSSQRT<TW) 

UW=<2.27*TW*SRTW/(TW+198.6) >*(10. )**(-8) 

SRGP=SQRT(2.*GBAR*( (PS-PINF) /RHOS) ) 

SRGK=SGIRT < 8H0W*UW*GB AR*SRGP/RFT ) 

Q=RENU*HS*U .-GO ) *SRGR/SIG 
OUDbsSRGP/RFT 
1N=0 
50 MN=0 
H=.U1 

NEQ=1 

NSO=U 

YY=0. 

X=0. 

DO *45 I=1»NP 

IF(l.NE.l) GO TO 32 

IND=4 

GO 10 34 

52 IND=1 

5b CALL RONKUT(NQEfNOS»H*XfFF»YY»YPPfFKKf AY Y , AYPP> SUMM1 » SUMM2 > IND ) 

59 NN=NM+1 

GM=l.-60GtNN) 

GM3=(GM)**3 

RH0G=1.-(GAMA1+GAMA2*GM3)*GM 

RH0U=SQRT( CRH0W*UW*GBAR*RFT)/(2.*SRGP) > 

FF= ( RHOG/RHOS ) *RHOU 

IND=IND+l 

IFUN.EQ.l) GO TO 39 

IFCNN.NE.IFP) GO TO 31 

DELTA=Y.Y 

IN=1 

GO TO 30 

5V IF(NN.NE.IG) GO TO 31 

DELTAT=YY 

ISDT=ISDT-1 

WRITE (b» lib) 

WR I TE ( 6 » 38 > DELTA » DELTAT » RENU » 0 * FLOG t H » I SOT 
58 FORMAT ( 4X.17H DELTA * 17H DELTA-T t 17H NU-RT 

T^RE JTTh 5 "7I7FT U(WS TT7R 5W 

.SIZE » 20H NO. ITERATIONS r/(2X»6E18.8rlOX. 13)) 

IF(NOFV.EQ.O) GO TO 99 
IFUCTR_.EQ.NOFV) GO TO 99 
GO TO 111 

31 1FUN0.LE.4) GO TO 3b 
4b CONUNUE 

STOP i 

> END 


NASA-Langley, 1968 12 S-121 
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